HARMONIC  MOTION  OF  CYLINDERS 
IN  UNIFORM  FLOW 


James  Thomas  Fry 


Monterey,  California  9394Q 


$  ffi 


PUS 


j  e  a 


onterey,  California 


in 


Hb 


*SSS«        VtEF' 


I    ■mi  ■—  — 


Harmonic    Motion    of 
Cylinders    in    Uniform    Flow 

by 

James  Thomas  Fry 


June  1975 


Thesi  s    Advi  sor 


T.    Sarpkaya 


:-"x  mnw  — 


Approved  for  public  release;  distribution  unlimited. 


T167976 


UNCLASSIFIED 


SECURITY   CLASSIFICATION   OF   THIS  PAGE   rWh.n  Date  Entered) 


REPORT  DOCUMENTATION  PAGE 


1.     REPORT   NUMBER 


2.  GOVT   ACCESSION   NO 


READ  INSTRUCTS. ns 
BEFORE  COMPLF.TI.Vf,  FORM 


3.      RECIPIENT'S  CATALOG   NUMBER 


4.     TITLE  (and  Subtitle) 


Harmonic  Motion  of  Cylinders 
in  Uniform  Flow 


5.     TYPE  OF   REPORT   ft   PERIOD   COVERED 

Master's  Thesis; 
June  1975 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORO; 

James  Thomas  Fry 


8.  CONTRACT  OR  GRANT  NUMBERf*; 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  CA   93940 


10.     PROGHAM   ELEMENT.  PROJECT     TASK- 
AREA  ft   WORK  UNIT  NUMBERS 


It.     CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  CA   93940 


12.     REPORT   DATE 

June    1975 


13.     NUMBER  OF  PAGES 


53 


14.     MONITORING  AGENCY  NAME  ft    AODRESSf//  different  from  Controlling  Office) 


Naval  Postgraduate  School 
Monterey,  CA   93940 


15.     SECURITY  CLASS,  (of  thle  t'oporl) 


UNCLASSIFIED 


15*.     DECL  ASSIFIC  ATI  OnTdOWN  GRADING 

SCHEDULE  S 


16.     DISTRIBUTION   ST  AT  EMEN  T  (of  thle  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.     DISTRIBUTION  STATEMENT  (of  the  abetract  entered  In  Block  20,  If  different  from  Report) 


18.     SUPPLEMENTARY  NOTES 


IS.     KEY  WORDS  (Continue  on  reveree  tide  It  necaeaary  and  Identify  by  block  number) 

Harmonic  motion 
Cable  st rummi ng 
Flow  about  oscillating  cylinders 


L 


20.     ABSTRACT  (Continue  on  reveree  aide  If  neceetary  and  Identity  by  block  number) 

The  time-dependent  force  acting  on  a  cylinder  undergoing 
harmonic  in-line  oscillations  in  an  otherwise  steady  flow  was 
measured  for  various  amplitudes  and  frequencies  of  oscillation 
and  mean  flow  velocity.   The  experiments  were  carried  out  in  a 
recirculating  water  tunnel  operating  as  an  open  channel  with  a 
free  surface  at  the  test  section.   The  time-dependent  force  has 
been  expressed  in  terms  of  a  mean  drag  coefficient  C,  and  the 


DD     i  JAN^S     1473  EDITION  OF    1  NOV  68  IS  OBSOLETE 

(Page      1)  S/N   0102-014-6601   ( 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Vhen  Dam  itntered) 


UNCLASSIFIED 


fliCUHITY    CLASSIFICATION    OF    THIS    PAGErH'.en    f>/-f«    £n(»f« 


20.   Abstract  (continued) 

Fourier-averaged  drag  and  inertia  coefficients  C.  and  C   as 
functions   of  the  relative  amplitude  A/D  and  the  frequency 
parameter  D/VT. 


DD     Form       1473 

,  1  Jan  73 
S/N    0102-014-CJG01 


UNCLASSIFIED 


SECURITY   CLASSIFICATION   OF   THIS  P»CEfl»fi«n  &*"   £"'■•"<*> 


Harmonic  Motion  of 
Cylinders  in  Uniform  Flow 


by 


James  Thomas  Fry 
Ensign,  United  States  Navy 
B.S.,  United  States  Naval  Academy,  1974 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
June  1975 


Library 

*  w»rforn/a  93940 


ABSTRACT 


The    time-dependent    force  acting     on    a    cylinder    under- 
going   harmonic    in-line    oscillations    in    an    otherwise    steady 
flow  was    measured    for    various    amplitudes    and    frequencies    of 
oscillation    and   mean    flow    velocity.      The   experiments    were 
carried   out    in    a    recirculating   water   tunnel    operating    as    an 
open   channel    with    a    free   surface    at    the    test   section.      The 
time-dependent   force    has    been    expressed    in    terms    of   a   mean 
drag    coefficient   C,    aand    the    Fourier-averaged    drag    and 

inertia    coefficients    C,    and    C      as    functions    of   the    relative 

d      m 

amplitude  A/D  and  the  frequency  paramter  D/VT. 
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NOMENCLATURE 


A  amplitude    of   cylinder    oscillation 

C  ■  mean    drag    coefficient 

C,  Fourier-averaged   drag    coefficient 

C  Fourier-averaged    inertia    coefficient 
m  3 

D  diameter    of   test    cylinder 

F  instantaneous    total    force    acting    on    the    test   cylinder 

Re  Reynolds    number    (Re    =    VD/v)   < 

T  period    of   oscillation 

t  time 

U  instantaneous    velocity 

U  maximum    velocity    in    a    cycle  (U      =    2-irA/T) 

m  J  J  K     m 

V  mean  velocity  of  the  uniform  flow 

e  phase  angle 

v  fluid    kinematic   viscosity 

p  fl ui  d    densi  ty 
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I.       INTRODUCTION 


Elastic    structures    of    one    or    more    degrees    of    freedom 
can   extract    energy    from    the    flow    about    them    and    can    develop 
catastrophic    flow-induced    vibrations.       The    understanding    of 
this    energy-extraction    process    is    of   paramount    importance 
if    one    is    either   to   eliminate    or   minimize    it    or    to    design 
the    elastic    structure    such    that    it    can   withstand    the    oscil- 
lations   under    the    contemplated   environmental    conditions. 

Even    a    superficial    familiarity   with    the    parameters 
involved    in    this    type    of    phenomenon    shows    that    the    investi- 
gation  of    the    flow-induced    vibrations    or   the    flow    about 
oscillating    bodies    is   more    than    an    extension    of    the    past 
studies    on    flow    about    bluff   bodies    at    rest.       In    fact,    the 
complexity    of   the    problem    is    increased   by    an    order   of 
magnitude.       In    view. of    this    fact,    it    is    no   wonder    that   the 
past    decade    has    produced   either   numerous    experimental    data 
with    the    goal    of    obtaining    ad-hoc    solutions    for   specific 
problems,   while    directing    little   if    any    attention    toward 
elucidation    of    the    underlying    mechanics    and    toward 
generalization    of   the    results    for    future    applications;    or 
numerous    analytical    models    most    or    all    of   which    had    nothing 
to    do   with    the    motion    of    the    very    medium   which    provided    the 
necessary    energy    to    set   the    body    in    motion.       In    fact,    the 
fluid   mechanics    of   the    flow-induced    oscillations    became    so 


incidental  to  the  phenomenon  that  the  models  dealt  essen- 
tially with  black-box-induced  oscillations.   To  be  sure, 
these  models  heavily  relied  on  experimental  data  partly  to 
justify  their  existence  and  partly  to  assign  numerical  values 
to  numerous  variable  constants  imbedded  in  them.   The  ability 
of  these  models  to  scale,  i.e.,  to  permit  extrapolation,  is 
uncertain.   Thus,  in  the  final  analysis  one  is  not  quite 
sure  whether  one  should  use  the  experimental  data  within  the 
range  of  their  application  or  the  curves  fitted  to  them  by 
the  empirical  models.   Suffice  it  to  say  that  the  phenomenon 
is  far  from  understood  even  for  the  idealized  conditions 
encountered  in  the  laboratory  without  the  alarming  conse- 
quences of  the  real  ocean  environment  where  fauna  and  flora, 
ever-changing  ocean  currents,  and  temperature  gradients  add 
further  complications  to  an  already  complex  problem.   Surely, 
over  design  is  not  the  answer  but  it  may,  in  the  next  decade 
or  so,  be  the  best  available  tradeoff  with  failure. 

Mathematical  models  of  flow-induced  vibrations  of  bluff 
bodies  and  the  response  of  circular  cylinders  to  vortex 
shedding  have  been  aptly  described  by  Parkinson  [1]  and 
Currie,  et  al  .  [2]  and  will  not  be  repeated  here. 

It  appears  that  among  the  various  models  considered  so 
far,  the  "wake-oscillator"  model  of  Hartlen  and  Currie  [3] 
attracted  more  attention  among  the  students  of  vibration 
analysis.   We  will  not  elaborate  here  on  the  attempts  made 
by  others  to  introduce  one  or  more  additional  terms  into  the 
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model  originally  proposed  by  Hartlen  and  Currie  for  such 
naive  attempts  produced  only  more  papers  and  unreal i sti cal ly 
defined  coefficients  than  sound  information  toward  the 
understanding  of  the  flow-induced  oscillations  both  in-line 
with  and  transverse  to  the  ambient  flow. 

The  in-line  oscillations  has  been  subjected  to  very 
little  investigation.   Chen  and  Ballengee  [4]  examined  the 
vortex  shedding  from  circular  cylinders  in  an  oscillating 
freestream  of  3  Hz  with  A/D  form  15  to  1,000,  D/VT  =  0.003, 
and  the  Reynolds  numbers  up  to  40,000,  the  vortex  shedding 
form  a  circular  cylinder  responds  instantaneously  to  the 
freestream  variations  and  that  "the  instantaneous  Strouhal 
number  stays  sensibly  constant  at  0.2+_0.01."   This  is 
rather  expected  since  the  amplitude  of  oscillations  is  many 
times  that  of  the  cylinder  diameter  and  the  flow  in  the 
absence  of  significant  accelerations  exhibits  a  quasi-steady 
behavior. 

Hatfield  and  Morkovin  [5]  studied  the  effect  of  an 
oscillating  freestream  on  the  unsteady  pressure  on  a 
circular  cylinder  (D/VT  from  0.15  to  0.25,  A/D  from  0.05  to 
0.087,  and  Re  =  50,000).   They  have  found  that  there  is  no 
significant  coupling  between  the  smal 1 -amp! i tude  freestream 
oscillations  and  the  vortex  shedding.   Their  results  would 
suggest  that  the  drag  coefficient  associated  with  the  mean 
flow  would  essentially  remain  constant  at  its  steady  state 
value.   This  study,  unlike  the  previous  ones,  concentrated 
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on  the  other  extreme  of  the  A/D  values.   In  this  case,  the 
flow  fluctuations  in  the  neighborhood  of  the  structure  are 
essentially  uncorre 1 ated  and  therefore  relatively  ineffec- 
tive.  Evidently,  for  values  of  D/VT  from  about  0.1  to  0.5 
and  A/D  from  about  0.2  to  1.0  that  the  interaction  of 
oscillations  with  the  body  becomes  significant  and  rather 
complex.   The  study  of  the  two  extremes  of  D/VT  and  A/D 
does  not  shed  much  light  on  the  understanding  of  oscillating 
flow  about  flow  about  bluff  bodies. 

Mercier  [6]  who  subjected  cylinders  to  large  streamwise 
oscillations  found  that  the  averaging  drag  coefficient 
significantly  increases  with  D/VT  and  that  the  rate  of 
increase  depends  on  the  amplitude  to  diameter  ratio  in  the 
range  of  A/D  from  0.2  to  3.0,  nD/V  from  0.1  to  infinity, 
and  Re  from  4,000  to  16,000. 

Davenport  [7]  subjected  bluff  bodies  (flat,  plate, 

circular  and  triangular  cylinders,  and  lattice  truss)  to 

small  amplitude  oscillations  in  a  water  flume  and  evaluated 

the  drag  and  inertia  coefficients  through  the  measurement 

of  the  rate  of  damping  of  the  amplitude  of  oscillations  of 

the  bodies,  i.e.,  without  directly  measuring  the  forces 

acting  on  the  bodies.   The  experiments  were  carried  out  in 

the  range  of  D/VT  from  0.02  to  2.0.   The  Reynolds  number 

ranged  from  about  1700  to  8,000.   Davenport  evaluated  C 

r  m 

and  C.  assuming  C,  to  be  equal  to  that  for  the  corresponding 
steady  flow  at  the  velocity  V  and  ignoring  the  term 
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involving  U  .   Since  C,  does  not  remain  constant  as  evidenced 
m         d 

2 
by  the  present  study  and  U   cannot  be  neglected  for  U  /V 

larger  than  about  0.1,  Davenport's  C.  and  C   values  are  not 

comparable  with  those  presented  herein. 

In  a  related  study,  Tseng  [8]  conducted  experiments 
with  flat  plates  normal  to  the  stream  undergoing  freely 
decaying  oscillations.   His  results  have  shown  that  the 
presence  of  the  mean  flow  significantly  increases  the  damping 
force  and  that  the  rate  of  extinction  of  the  oscillation 
increases  monotoni cal ly  with  speed. 

Tanida,  Okajima,  and  Watanabe  [9]  found  that  for  a 
circular  cylinder  oscillating  parallel  to  the  flow 
(A/D  =  0.14,  D/VT  from  zero  to  0.5,  and  for  Re  =  80  and 
4,000),  the  vortex  synchronization  can  be  observed  in  a 
range  around  double  the  Strouhal  frequency,  where  vortices 
are  shed  with  a  frequency  half  the  imposed  one.   Tanida, 
et  alli's  results  show  that  the  mean  drag  reaches  its  maximum 
in  the  middle  of  the  synchronization  range,  i.e.,  D/VT 
between  0.2  and  0.4,  and  that  the  sign  of  C.  is  such  that  no 
energy  can  be  extracted  from  the  fluid  to  render  the  oscilla- 
tions unstable.   In  other  words,  in-line  oscillation  is 
stable  for  the  two  Reynolds  numbers  tested.   Tanida  conjec- 
tured that  instability  is  likely  to  occur  at  much  higher 
Reynolds  numbers.   Unfortunately,  their  C.  values  cannot  be 
relied  upon  since  the  inertial  forces  were  subtracted  only 
approximately  and,  according  to  the  data  reported  herein, 
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are  far  from  being  correct.   Thus,  their  C,  values  cannot  be 
considered  sufficiently  accurate  to  assess  the  stability  of 
i  n-1 i ne  osci 1 1 ati  ons  . 

Goddard  [10]  carried  out  a  numerical  solution  of  the 
drag  response  of  a  cylinder  to  streamwise  velocity  fluctua- 
tions for  Re  =  40  and  nD/V  =  0.019,  0.12,  and  3.18,  and  for 
Re  =  200  and  nD/V  =  0.149.   He  found  that  for  very    low 
frequencies  the  instantaneous  values  of  drag  correspond  very 
nearly  to  the  quasi-steady  solution  and  that  for  higher 
frequencies  the  drag  anticipates  the  freestream  velocity 
maximum.   This  work  cannot  be  generalized  to  higher  Reynolds 
numbers  since  it  is  based  on  the  Navier-Stokes  equations  and 
since  the  diffusion  of  vorticity  in  the  concentrated  vortices 
for  Re  larger  than  about  200  is  primarily  turbulent. 
Numerical  or  exact  solutions  of  the  Navier-Stokes  equations 
cannot  take  into  consideration  such  a  turbulent  diffusion 
unless  modified  through  the  use  of  an  appropriate  eddy 
v  i  s  c  o  s  i  ty  . 

DECOMPOSITION  OF  MEASURED  FORCES 

The  numerical  calculations  as  well  as  measurements  in 
time-dependent  flow  yield  the  resultant  force  as  a  function 
of  time  for  a  given  set  of  numerical  values  of  the  independent 
parameters.   Thus  it  is  not  possible  without  a  suitable 
hypothesis,  to  express  the  force  both  as  a  function  of  time 
and  remaining  independent  parameters  as  one  ordinarily  would 
in  a  closed  form  solution.   Such  a  working  hypothesis  is 
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particularly  necessary  for  the  cable  strumming  problem  since 
the  results  are  to  be  incorporated  into  the  dynamics  of  the 
cylinder  motion  in  the  form  of  a  forcing  function.   It  should 
be  stated  at  the  outset  that  there  is,  at  present,  no 
generally  accepted  hypothesis  to  decompose  the  time-dependent 
force  into  suitable  components.   As  it  will  be  seen  shortly, 
even  the  existing  hypotheses,  such  as  the  so-called  Morison's 
equations,  are  not  applicable  to  periodic  flows  with  a 
non-zero  mean  velocity. 

Stokes,  in  a  remarkable  paper  on  the  motion  of  pendu- 
lums, showed  that  the  expression  for  the  force  on  a  sphere 
oscillation  in  an  unlimited  viscous  fluid  consists  of  two 
terms,  one_ involving  the  acceleration  of  the  sphere  and  the 
other  the  velocity.   This  analysis  shows  that  the  inertia 
coefficient  is  modified  because  of  viscosity  and  is  augmented 
over  the  theoretical  value  valid  for  irrotational  flow.   The 
drag  coefficient  associated  with  the  velocity  is  modified 
because  of  acceleration,  and  its  value  is  greater  than  it 
would  be  if  the  sphere  were  moving  with  a  constant  velocity. 
In  general,  the  force  experienced  by  a  bluff  body  at  a  given 
time  depends  on  the  entire  history  of  its  acceleration  as 
well  as  the  instantaneous  velocity  and  acceleration.   Thus, 
the  drag  coefficient  in  unsteady  flow  is  not  equal  to  that 
at  the  same  instantaneous  velocity  in  steady  flow.   Neither 
is  the  inertia  coefficient  equal  to  that  found  for  unseparated 
potential  flow.   As  yet  a  theoretical  analysis  of  the  problem 
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for  separated  flow  is  difficult  and  much  of  the  desired 
information  must  be  obtained  both  experimentally  and 
numerically.   In  this  respect,  the  experimental  studies  of 
Morison  and  his  co-workers  [11]  on  the  forces  on  piles  due 
to  the  action  of  progressive  waves  have  shed  considerable 
light  on  the  problem.   The  forces  are  divided  into  two  parts 
one  due  to  the  drag,  as  in  the  case  of  flow  of  constant 
velocity,  and  the  other  due  to  the  acceleration  or  decelera- 
tion of  the  fluid.   This  concept  necessitates  the  introduc- 
tion of  a  drag  coefficient  C.  and  an  inertia  coefficient  C 
3  d  m 

in  the  expression  for  force.   In  particular  if  F  is  the 
force  per  unit  length  experienced  by  a  cylinder,  then 


F ■=  0.5  C.pDIUIU  +  C  pirD^/4  dU/dt 
d  m 


where  U  and  dU/dt  represent  respectively  the  undisturbed 

velocity  and  the  acceleration  of  the  fluid. 

On  the  basis  of  irrotational  flow  around  the  cylinder, 

C   should  be  equal  to  2  (cylinder  at  rest,  the  fluid 

accelerating;  otherwise  C   =  1),  and  one  may  suppose  that 

the  value  of  C,  should  be  identical  with  that  applicable  to 

a  constant  velocity.   However,  numerous  experiments  show  that 

this  is  not  the  case  and  that  C,  and  C   show  considerable 

d  m 

variations    from    those    just    cited    above.      Even    though    no    one 
has    suggested    a    better    alternative,    the    use    of   the    Morison's 
equation    gave    rise    to    a    great    deal    of    discussion    on   what 
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values  of  the  two  coefficients  should  be  used.   Furthermore, 
the  importance  of  the  viscosity  effect  has  remained  in  doubt 
since  the  experimental  evidence  published  over  the  said 
period  has  been  quite  inconclusive. 

The  drag  and  inertia  coefficients  obtained  from  a  large 
number  of  field  tests,  as  compiled  by  Wiegel  [12],  show 
extensive  scatter  whether  they  are  plotted  as  a  function  of 
the  Reynolds  number  or  the  so-called  period  parameter  U  T/D. 

The  reasons  for  the  observed  scatter  of  the  coefficients  C 

m 

and  C,  remained  largely  unknown.   The  scatter  was  attributed 
to  several  reasons  or  combinations  thereof  such  as  the 
irregularity  of  the  ocean  waves,  free-surface  effects, 
inadequacy  of  the  average  resistance  coefficients  to  repre- 
sent the  actual  variation  of  the  nonlinear  force,  omission 
of  some  other  important  parameter  which  has  not  been 
incorporated  into  the  analysis,  the  effect  of  ocean  currents 
on  separation,  vortex  formation,  and  hence  on  the  forces 
acting  on  the  cylinders,  etc. 

The  most  systematic  evaluation  of  the  Fourier-averaged 
drag  and  inertia  coefficients  has  been  made  by  Keulegan  and 
Carpenter  [3]  through  measurements  on  submerged  horizontal 
cylinders  and  plates  in  the  node  of  a  standing  wave,  applying 
theoretically  derived  values  for  velocities  and  accelerations 
Additional  measurements  have  been  made  by  Sarpkaya  [14]  of 
the  in-lines  as  well  as  transverse  forces  acting  on  cylinders 
and  spheres  in  a  sinusoidal ly  oscillating  fluid  and  it  was 
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found  that  the  drag  coefficient  as  well  as  the  inertia 

coefficient  for  a  strictly  sinusoidally  oscillating  fluid 

(no  mean  velocity)  is  a  function  of  U  T/D  and  that  the 

effect  of  the  Reynolds  number  is  rather  secondary  and 

certainly  obscured  by  the  excellent  correlation  of  the  data 

with  the  period  parameter  U  T/D. 

On  the  basis  of  the  above  discussion,  one  would  assume 

that  Morison's  equation  would  apply  equally  well  to  periodic 

flow  with  a  mean  velocity  where  u  =  V  -  U  cose  and  that  C, 

m  d 

and  C   will  have  constant,  time -invariant,  Fourier  or  least- 
m 

squares  averages.   This,  inturn,  implies  that  C,  and  C   are 

independent  of  the  associated  flow  phenomena.   There  is, 

however,  no  a  priori  assurance  in  the  principles  of  fluid 

mechanics  of  theory  of  models  that  this  is,  in  fact,  the 

case.   Thus  the  effect  of  the  combination  of  a  uniform 

current  and  harmonic  oscillations  on  the  time-average  and 

oscillatory  forces  acting  on  circular  cylinders  will  have 

to  be  re-examined  and  the  limits  of  application  of  the 

Morison's  equation  be  delineated. 

It  is  a  priori  evident  that  both  u  =  -U  cose  and 
r  m 

u  =  V  -  U  cose  yield  the  same  acceleration  du/dt.   Thus, 
the  force  in-phase  with  the  acceleration  in  Morison's 
equation  remains  unaffected  by  the  presence  of  the  mean  flow 
The  results  presented  herein  show  that  this  is  not  the  case. 
Furthermore,  the  use  of  the  Morison  equation  as 
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F   =   0.5pC.(tf-U    coso)  ]  V-U    cose  I    +   C   PttD2/4.    du/dt 
dm  '         m  '  m  ' 

requires    that    the    time-averaged    drag    force    be    calculated    by 
increasing    the    force   calculated    from    the    steady    flow   by    a 
factor    [1    +   0.5(U    /V)    ].      The    results    presented   herein    show 
that    such    an    analysis    appreciably    underestimates    the    measured 
mean    forces.      It    suffices    to   state    that   the    fluid    flow 
phenomena    for    bluff   bodies    are    significantly    affected    by    the 
combination    of   currents    and    harmonic    oscillations    and    that 
the    results    for   steady    currents    alone    and   oscillations    alone 
cannot    be    combined   to  yield    reliable    estimates    of   forces    due 
to   both    acting    together. 

The    time-dependent   forces    in    the    present   study    are 
analyzed    according    to   the    following    three-coefficient 
equati  on 


F   -    0.5CdpDV2    +   Cmp^^(-Umcos^t) 


2tt 


27f 


CdDUmp£"  I  cosTLt  jcos^pt 


which   may   be   written    as 


0.5PDV 


12    =    Cd    +    CmTT2(UmT/D)(D/VT)2sin^t    - 


2tt 


Cd(UmT/Dr(D/VTr    |cos^-t|    cos^-t 


19 


in  which  C   and  C,  are  given  by  their  Fourier  averages  as 
m      d 


2tt 


Cm    =    (2UmT/7T3D)  J~    (F   sine)    de/(PU^D) 


and 


2tt 
Cd   =   -(3/4)  J"    (F   sine)    de/(PU^D) 


Evidently,  C.  ,  C  .  and  C,  are  functions  of  VT/D  and 
U  T/D  or  A/D.   They  may  depend  also  on  the  Reynolds  number 
which  does  not  explicitly  appear  in  the  above  expression 
because  of  the  assumptions  made  in  the  formulation  of  the 
basic  force  equation. 

In  the  foregoing,  neither  the  coefficient  C,  is  assumed 
to  be  equal  to  the  steady-state  drag  coefficient  for  a 

uniform  flow  at  the  constant  velocity  V,  nor  C   and  C,  are 

J  m      d 

assumed  to  be  identical  to  those  obtained  for  a  strictly 

harmonic  oscillation.   In  fact,  the  results  show  that  C,  = 

C,  (steady)  only  for  U   =  0,  and  C,  and  C   are  equal  to 
d       j  i  j  m  d      m      ^ 

those  obtained  for  the  harmonic  oscillation  only  for  VT/D  =  0 
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II.   EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


A.    EQUIPMENT 

1 .    NPS  Water  Tunnel 


The  experiments  were  performed  in  a  recirculating 
water  tunnel  (Figure  1)  which  had  a  capacity  of  approximately 
500  gallons.   The  galvanized  test  section  was  four  inches 
wide,  eight  inches  high,  and  sixteen  inches  long.   A  low-rpm, 
high-capacity,  fourteen-inch-di ameter-di scharge  centrifugal 
pump  was  used  to  circulate  the  fluid  through  the  test 
section.   The  velocity  of  the  fluid  was  regulated  by  a 
butterfly  valve  arrangement  which  was  situated  downstream  of 
the  test  specimen.   Velocities  of  .7  to  1.5  fps  were  obtained 
by  adjusting  the  vains  of  the  butterfly  valve. 
2  .    Pendulum  and  Motor  Arrangement 

The  harmonic  motion  in  this  experiment  was  obtained 
by  a  small,  variable  speed,  electric  motor  and  crank  arrange- 
ment (Figure  2).   The  crank  was  attached  to  the  arm  of  a 
pendulum  41  inches  from  the  pivot  point.   The  overall  length 
of  the  pendulum  was  56  inches.   The  frequency  of  the  pendulum 
oscillation  was  regulated  by  the  speed  of  the  electric  motor. 
The  frequency  of  oscillation  varied  from  0.6  to  9.0  c.p.s. 
The  amplitude  of  the  oscillation  was  set  by  adjusting  the  end 
of  the  crank  on  a  rotating  disk  (Figure  2).   The  amplitude 
"A",  ranged  from  3/16  to  1.5  inches. 
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3 .    Test  Specimen 

Circular  cylinders  of  3/4  inch  and  one  inch  were 
used  as  the  test  specimen.   The  3/4  inch  cylinder  was  made 
of  plexiglass  and  the  one  inch  cylinder  was  made  of  aluminum 
tubing.   The  test  specimen  was  switched  to  aluminum  tubing 
because  the  inertia!  force  experienced  was  less  than  that 
experienced  with  the  plexiglass  cylinder. 

The  cylinders  were  attached  to  the  arm  of  the 
pendulum  by  means  of  an  aluminum  yoke  arrangement  (Figure  3). 
The  action  of  the  pendulum  caused  the  cylinder  to  oscillate 
in  the  center  of  the  test  section  and  in-line  with  the 
stream  flow. 

In  order  to  oscillate  the  cylinder  with  the 
pendulum  arrangement,  it  was  necessary  to  operate  the  tunnel 
with  the  top  of  the  test  section  open.   The  water  level  was 
maintained  as  close  to  the  top  of  the  test  section  as  possible 
in  order  to  make  the  free  surface  effects  negligible. 
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STRAIN 
GAG-ES 


SCALE-  ONE    HALF  DIMENSIONS:    INCHES 

Figure  3.   Aluminum  yoke/force  transducer 
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D    =    1 .0    INCHES 

A/D    =    0.5 

T    =    0.485     sec 


IN-LINE    FORCE 


WET 


DRY 


ACCELERATIOI 


Figure  4.   Force  and  acceleration  traces. 
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D    =    1  .0    INCHES 

A/D    =    0.5 

T    =    0.162    sec 


IN-LINE    FORCE 


WET 


DRY 


ACCELERATION 


Figure  5.   Force  and  acceleration  traces. 
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4 .    Sensor  and  Sensor  Support 

The  mean  velocity  of  the  fluid  was  derived  from 
the  pressure  differential  of  a  pitot  tube  located  upstream 
of  the  test  specimen. 

The  instantaneous  acceleration  of  the  pendulum  arm 
was  monitered  by  an  accelerometer  attached  to  the  arm  of  the 
pendu 1  urn . 

The  in-line  force-measuring  device  consisted  of  a 
cantilever  beam  as  part  of  the  yoke  arrangement.   Four 
piezoresi s ti ve  strain  gages  were  mounted  on  the  cantilever 
beam  and  properly  waterproofed.   A  drawing  of  this  force 
transducer  is  shown  in  Figure  3.   The  transducer  was 
calibrated  by  hanging  loads  at  the  mid  section  of  the 
cylinder  in  the  vertical  direction. 

B.    PROCEDURE 

It  was  realized,  during  the  early  stages  of  the 
experiment,  that  a  procedure  of  recording  both  the  wet  and 
dry  oscillations  at  a  common  frequency  would  be  necessary 
for  a  complete  evaluation  of  the  data.   The  wet  oscillating 
force  was  made  up  of  an  inertial  force  and  a  resistive  force 
produced  by  the  fluid.   In  order  to  separate  the  forces,  the 
dry  force  was  recorded  at  the  same  frequency  as  the  wet 
force  and  then  subtracted  from  the  wet  force  to  give  the  pure 
resistive  force. 
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To  accomplish  the  physical  recording  of  the  wet  and 
dry  force  oscillations,  the  water  level  in  the  tunnel  was 
raised  and  lowered  for  each  frequency  tested.   The  proper 
status  of  the  cylinder  was  noted  on  the  chart  recordings. 

Figures  4  and  5  show  the  wet  and  dry  forces  as  recorded 
for  two  typical  frequencies.   The  acceleration  curve  was 
used  as  a  reference  in  matching  the  wet  and  dry  run  curves. 
It  was  convenient  that  the  acceleration  traces  were  nearly 
perfect  sine  curves. 

From  traces  similar  to  the  one  shown  in  Figure  4,  the 
in-line  force  was  read  and  recorded  from  both  the  wet  and 
dry  curves.   The  result  was  then  punched  on  IBM  computer 
cards  for  every    1.0  or  2.0  mm ' s  . 

Finally  the  drag  and  inertia  coefficients  were  calcula- 
ted through  the  computer  program  given  in  Appendix  A. 
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III.   DISCUSSION  OF  RESULTS 

The  results  of  the  present  investigation  will  be 
reported  in  terms  of  C.,  C  ,  ,  and  C   as  a  function  of  D/VT 
or  VT/D  for  various  values  of  A/D.   Figures  6  through  9  show 
C.  as  a  function  of  D/VT  for  A/D  =  0.25,  0.5,  1.0,  and  1.5. 
In  all  cases  C,  is  equal  to  its  steady-state  value  for 
D/VT  =  0,  i.e,  for  f  =  0.   The  mean-drag  coefficient  reaches 
a  maximum  value  at  D/VT  between  0.2  and  0.35.   The  value 
of  D/VT  at  which  this  maximum  occurs  shifts  from  0.35  to  0.2 
as  A/D  increases.   Furthermore,  the  maximum  value  of  C, 
increases  with  increasing  A/D.   For  example,  for  A/D  =  1.5 
and  D/VT  =  0.1  (a  wave  of  15  ft  amplitude  and  20  second 
period  superimposed  on  a  current  of  5  ft/sec  velocity  acting 
on  a  10  ft  diameter  cylinder),  the  mean  drag  acting  on  a 
cylinder  is  increased  by  nearly  50% 

Upon  reaching  a  maximum,  C,  decreases  sharply  and  as 
evidenced  by  the  results  of  A/D  =  0.25,  reaches  a  value 
nearly  equal  to  that  corresponding  to  a  steady  flow  at  the 
velocity  V.   Thus,  for  a  given  set  of  V  and  D,  in-line 
oscillations  of  either  very    small  or  very    large  frequencies 
do  not  alter  the  mean  resistance,  e.g.,  for  A/D  =  0.25,  C. 
remains  essentially  constant  in  the  range  0.45<D/VT<0.1  or 
in  the  range  2>VT/D>10.   In  other  words,  when  the  frequency 
1/T  =  f  exceeds  twice  the  Strouhal  frequency,  there  is  no 
more  vortex  synchronization. 
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From  a  practical  point  of  view,  the  increase  of  C\  to 
values  several  times  larger  than  the  steady-state  value  is 
quite  significant.   Field  studies  conducted  to  determine 
wave  forces  on  cylinders  may  be  strongly  affected  by  the 
presence  of  mean  currents  and  may  yield  drag  and  inertia 
coefficients  which  are  significantly  different  from  those 
obtained  with  purely  harmonic  oscillations.   It  is,  of 
course,  fully  realized  that  the  ocean  waves  are  not  usually 
simple  harmonic,  but  are  rather  irregular. 

Figures  6  through  9  also  show  two  additional  facts 

of  special  importance.   Firstly,  as  far  as  C,  is  concerned, 

there  is  no  special  significance  of  the  case  where  U  >V.   A 
r  m 

simple  calculation  shows  that  U  >V  for  D/VT>0.637  for 

m 

A/D  =  0.25;  for  D/VT>0.318  for  A/D  =  0.5;  for  D/VT>0.159 

for  A/D  =  1.0;  and  D/VT>0.106  for  A/D  =  1.5.   Secondly,  C, 

is  not  equal  to  Cd  (steady)  ( 1 +0 . 5U2/V"2) .   In  fact,  only 

for  wery    small  values  of  U  /V  that  the  time-averaged  drag 

m  33 

force   can    be    calculated   by    increasing    the    force    calculated 

from    the    steady    flow    by    a    factor    ( 1 +0 . 5U2/V~2)  . 

m 

•     The    Fourier-averaged    drag    and    inertia   coefficients    C. 
and    C      are    shown    in    Figures    10    through    13    and    14    through    17 

m  3  3  3 

respectively  as  function  of  VT/D  for  A/D  =  0.25,  0.5,  1.0, 

and  1.5.   The  coefficients  C,  and  C   have  also  been  calculated 

d  m 

through    the    use    of   the    least-squares    method    and    a   modified 
least-squares    method.      Since    the    results    differ      wery    little 
from   those    obtained    through    the    use    of   the    Fourier   analysis, 
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only  the  Fourier-averaged  coefficients  are  presented  herein. 

On  each  graph,  the  values  corresponding  to  VT/D  =  0  are 

taken  from  Sarpkaya  [14]. 

Firstly,  it  is  noted  that  both  C,  and  C   exhibit  larger 

dm 

scatter    than    C..      The    reasons    for    this    scatter    are    rather 
d 

understandable.       For    a    given    A,    D,    and    V,    A/D    is    fixed    and 
only    T    is    variable    in    VT/D.      Thus,    yery    small    values    of 
VT/D    correspond    to    yery    small    values    of   T    or   to    rather    high 
frequencies    of    oscillation.      This    means    that   not   only    the 
fluid    forces    but    also    the    inertia!    forces    due    to   the   mass    of 
the    oscillating    system    are    large.      Furthermore,    at    very    high 
frequencies,    the    acceleration    is    not   exactly    harmonic    and 
the    second    order    accelerations    may    introduce    some    error    into 
the    evaluation   of    the    data.      In    passing,    it   should   be    noted 
that    attempts    to    increase    T    by    decreasing    V    so    as    obtain    the 
same    value    of    VT/D    are    not    very   desirable    since    this    leads 
to   smaller    Reynolds    numbers.      This    and   similar   other    facts 
simply    point    out    the    extreme    difficulty    of    experimentation 
in    unsteady    flows    and    are    sufficient    to    put    the    critical 
reader    in    touch   with    his   wisdom. 

For    large    values    of    VT/D,    the    frequencies    are    very    low 
and   give    rise    to    rather   small    fluctuating    forces    about    the 
mean    drag    force.      Thus,    the    decomposition    of   the    small 
alternating    component    of   the    total    force    into    in-phase   and 
out-of-phase    components,    i.e.,    into   drag    and    inertial    forces, 
leads    to    a    number    of    difficulties.       In    this    range,    small 
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errors  in  phase  between,  say,  zero  acceleration  and  the 
actual  position  of  the  cylinder  in  the  cycle  lead  to 

relatively  large  errors  in  the  evaluation  of  C,  and  C  .   In 

J  3  dm 

spite    of   this,    however,    the    data    are    sufficiently    consistent 
to    draw   quantitative    conclusions.      It    is    apparent    from 
Figures    10    through    13    that    C.    increases    rapidly   with 
increasing    values    of   VT/D    and    remains    positive    throughout 
most    of   the    VT/D    values.       Only    for    A/D    -   0.25    and    0.50    and 
for   \jery    small    values    of   VT/D    that   C.    becomes    negative. 
This    range    of    VT/D   values    nearly    fall    in    the    velocity    range 
where    U    >V.      Evidently,    the    cylinder   cannot   extract   energy 
from   the    fluid   when   C  ,>0  .      Thus,    the    negative    values    of   C, 
in    a   small    region   of   VT/D    values    imply    that    it    is    possible 
to   excite    in-line    oscillations    of   very   small    amplitudes    at 
frequencies    five    or   more    times    the    Strouhal    frequency. 
Such    an    oscillation   has    not    previously   been    observed    partly 
because    no    one    has    excited    the    elastic   test    cylinders    at 
such    high    frequencies    and   partly   because    the    cylinders    were 
allowed   to    oscillate    in    all    directions    with    the    consequence 
that    in-line   oscillations    may   have    been    obscured   by    rela- 
tively   larger    transverse    oscillations. 

The    inertia    coefficient   C      is    shown    in    Figures    14 

m  3 

through  17  as  a  function  of  VT/D  for  various  values  of  A/D. 
For  small  values  of  VT/D,  C   approaches  its  ideal  potential- 

m    rr  r 


flow  value  of  unity.   Near  Strouhal  frequency,  however,  C 
decreases  sharply,  particularly  for  small  amplitudes  of 


m 
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oscillation.   In  other  words,  near  Strounal  frequency,  the 
force  acting  on  a  cylinder  undergoing  small  amplitude 
in-line  oscillations  is  primarily  from  drag.   This  is  rather 
interesting  if  we  consider  the  fact  that  for  small  values 
of  A/D,  say  A/D  =  0.25  or  UmT/D  =  1.57,  the  force  acting  on 
a  cylinder  oscillating  in  a  fluid  otherwise  at  rest  (i.e., 
V  -  0)  is  essentially  of  inertial  nature  [14].   Thus,  the 
presence  of  a  mean  flow  superimposed  on  harmonic  oscilla- 
tions can  significantly  alter  not  only  the  magnitude  of 
the  inertial  force  but  also  the  nature  of  the  entire  flow 
structure  about  the  body.   For  large  values  of  VT/D  where 

U/V  is  rather  small  (U  /V  =  0.15  for  A/D  =  0.25,  and 
m  v  m 

U  /V  =  0.30  for  A/D  =  0.50),  the  forces  acting  on  the  body 
m 

follow  the  changes  in  velocity,  as  noted  in  connection  with 
the  discussion  of  C  , ,  and  the  inertial  component  of  the 
force  is  negligibly  small.   However,  it  must  be  emphasized 
once  again  that  the  quasi-steady  model  is  grossly  inadequate 
to  predict  the  average  drag  and  the  response  of  the  instan- 
taneous force  to  changes  in  velocity  for  large  periods  or 
for  extremely  small  frequencies  bears  no  relation  to  the 
existence  of  an  actual  quasi-steady  state  of  flow. 
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IV.   CONCLUSIONS 

The   experimental    investigation    of   the    in-line    oscilla- 
tions   of    a    circular   cylinder   in    a    flow  with    a   mean    velocity 

V   has   yielded   the    force    coefficients    C,,    C    ,    and   C,    and    has 
J  dm  d 

shown    that: 

(a)  The  mean  flow  has  significant  effects  on  C,,  C 

v    '  3  dm 

and  C,  and  that  the  results  of  the  experiments  with  harmonic 
oscillations  in  a  fluid  otherwise  at  rest  are  not  applicable 
to    oscillations    of   a   cylinder    in    uniform    flow; 

(b)  In    a    region    of   small    values    of   D/VT  where    the 
frequency    of    oscillations    is    five    or   more    times    larger   than 
the    Strouhal    frequency,    C,    is    negative    and   energy   may   be 
transferred    from    the    fluid   to   the    cylinder.      This,    in    turn, 
may   give    rise    to   self-excited    in-line    oscillations    in    a 
manner   similar   to   those    for   the    transverse    oscillations;    and 
that 

(  c)      The    most    significant    region    of   the    D/VT   values    in 
which    the    force    coefficient   differ   significantly    from    those 
obtained    in    steady    flow    or   with    harmonic    oscillations    in    a 
fluid    at    rest    occur    in    the    vicinity    of    the    Strouhal    numbers. 
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APPENDIX  A 
COMPUTER  PROGRAMS 


cc 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

cc 


cccccccccccccccccccccccccccccccccccccccccccccccc 

AVERAGE    CD    AMD    CM    CALCUL ATir.NS    f 'OP    CYLINDERS 
TIPL=OI?lEh?!CNLLS    TIME     (  T I  ME/PE"  IOD  J 

bet;  =ui  iEr;s ic;r  li  s  di  splacement  (umax*per/di  a) 

CL=CYLI'Ji)Ec<    LENGTH    IN    FL"FT 
DIA=D!  AM.ETEP    OF    CYLIM0E0     IN    FEET 
AMP=AMPL  ITUOE     Df     MfTiorj     ifj    rEET 
PER-PERITD    IN    SCCS    (CHART    PERIOD/Cl 
ZI=FCRCE     COEFFICIENTS 


'APTS"EED) 


CM=INEP^!A  COEFFICIENT 
CO = DRAG  COEFFICIENT 
REMF=REMAINOER  FUNCT  ION 
A  1 1  B  1=  Ff  UP  I  ER  CO  E  F  F  t  C  I  ENTS 
N=NUMBER  OF  DA"^  SF"rs 
CN=CCNVERTIC.N  FACTOR 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


DIMENSION  FCPCE( 100) 
READ  IN  PAPAMETEPS  WHICH  CHARACTERIZE  THE 


DAT 


CCCCCCCCC 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
cccccc:cc 

A  SET 


-   N  =  12 

DO    100    I*1,N 
RE  AD  (  5 , 1 0  J  DI  A  i  AN"  P  ,  P E  P. ,  NC  fi  P  D ,  UM X  ,  CM  t  F MMM ,  NRUN 

READ (5 t 55)     (FCRCE(K),     K=1,NCAP.D) 

--G-32.17^  -  - 

RH0=62.4/G 
PI=3.1'tl59 

cnu=c.  oc  00105 

CL^0.33:J?334 
_  V3AF  =0,9563 

INITIATE    DATA    SETS 
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85 


CCMP 

Zl=( 

Z2=( 

Z2LS 

CCMP 

BET.' 

P.EYN 

WRIT 

WRIT 

DO    t 

SUM 

SUM 

CCN1 

CO  PR 

DO    8 

FOPC 

FOPC 

c:n- 

FMNM 
FMVM 
TIME 
CM=0 
CD  =  0 
CMLS 
CDLS 
FAA- 
FKB  = 
FCC^ 
VLO^ 
F  E  E  = 
DELT 
J  =  NC 
DO  2 
F  =  C\ 
ALPH 
SINA 


FORCE    CC  EFFICIENTS 
R*  PER  )  /{  P I  *P  I*PT  *0I  A*D  I  A-C L  *P.HO* 
ER*PER)/  (  3-r  HC*D  !  A*°  I*CL*AMP*AM'"> 
*PEF*PFt-  )  /  (3*PI*P!*RHC*,)IA*CL*AM 

RF~A 

I  *  A  r ■'  P  /  D T  A 

BAR*DIA)/CNU 

15  )0i  A  ,.-  MP,  PEPiUNXiCNf  \'RUN 

20) 


UTE 

^  *  •  f 

=  (-  + 

UTE 
=  "> #  p 

0=(V 
E(6, 
E(6, 

5  K= 1,  NCAPD 

=  0.0 
=  S  U  -1 
INUE 


AMP) 

) 

P*AMP) 


F  (  K  } 

E  {  K  ) 

INUE 

-  f 

=  0.0 
.0 
.0 
=  0 
=  0 
0.0 
0.0 
0.0 

o.o 

0.0 
AT  •-  L 

AP.D- 

0  0  K 
*  (  F  " 
A=«>* 

=  S  I M 


+  FCRCF (K) 

UM /NCAPD 
1,  NCARD 

=  FCPCE(K)  -  CORP. 

=  -  FORCE (K ) 


MM  M 
FMMM 


-  CPRR 


.O/NCARD 

1 

=  itJ 

PCE(K)  +1  f  RCE  (K  +  l  )  )/2.0 

Pl*<  "I  Mi  H)[[  TAT/2.0) 

(ALPHA) 
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cn$A*cos 

f  SIMA  = 

fci,:a= 
fls^ 

I  LS  = 

F/ AS=2*P 
Fri  155=2*" 
F( CS  =  2* 
F0DS=2 


(  Al  PHA) 

F*f  I  MA 

F*CFSA 

r*Ct)SA*AUS< 

f*co:  A*/.rtS( 
I  "F*l  ,onoo*CPS 
I  *f  *f  *Ci   SA*ABS 

'M*f>"l !':;- s:r  a* 


Cn 

SA1 

CO 

LA) 

A- 

rnsA 

*cr 

SA*C 

'"SA*^ 

(C 

r-  S  A  ) 

*  l . 

f  00 

*F 

AP 

S(C'.- 

SA) 

n.o 

CCO*F 

L. 

A0  1  j 

-F 

00 

0*F 

tCAPD)  1/2. 
VP-.O) 


FEES  =2*'>I  *F*F*S!NA*1  .0 
FAA=FAAS+FAA 
FBl3=F3BS+FBB 
FCC=Ff  cs+rcc 
FOC=FUPS-M  00 
FEE=FEES+FEF 
C,"4=FSIN\+CM 
Cf)  =  f  COSA+CO 
CMLS=FSINA    +CMLS 
CDLS=FLS+COLS 

V»RI7E(  6,  3  J)'.'  I  ME,  ALPM.A,COSAfSINA|Ff  FCCSA,  FSIMA 
T1HF=T!  4F«-DLLTA~ 
200    CCJNTINUc 

T!f'F=TI  1F+DELTAT 
F  =  CSi*(FLiRCE  (  1  )+r^PCE(! 
Al  P  IA=2-r>!  -J  (TT,VL  +  Dti  T 
SINA=SIN(ALPHA) 
COSAsCCS ( ALPHA) 
FSThA-  F*SINA 

FCCSA=  F*CCSA 

FAAS=2*P  !*r  -:1  .0'-  GO-COS 
F  p  f>  s  =  2  *  P  I  *  F  *  F  *  0  S  A  *  /  B  S 
FCCS=2*°  !*r*5TMA*CfJSA* 
FDDS=2*P!*F*S IW  *S IMA* 
FEES=2*PI*F*f  *SIf:A*l  .0' 
FAA=FAAS+FAA 

FBB=FBBS+FBB  . 

FCC=FCCS+FCC 
FDC=FODS+FOD 
FEt-FFES+FEE 
CM=FSINV+CM 

co=rcnsA+co 

CDLS=FLS+CDt S 

WR I TE ( 6, 30 ) T I  ME, ALPHA , CCSA , S I N A , F , FCCS A , FS I M  A 

CM=Z1*CM  •-:=DrLTA~ 

CD=Z2-C0  *0£ LTAT 

CMLS=Zl*Cf  LS    *PFLTAT 

CDl  S-=Z2LS*CI)LS    *OELTAT 

CMFF=( Z1/2.0)*(FEE*FAA-FCC*F3B) /(FCD-FA/ -FCC -FCC ) 


A* 

COS  A 

*cn 

SA* 

ens  a 

(C 

OS  A) 

*i. 

000 

*F 

AB 

S(CH 

SA  ) 

*1. 

0000 

1. 

oooo 

0*F 

*F 

0', 

CDFT 

WR! 

WR!T 

/  NGl 

W  IT 

TIME 

DO    3 

F  =  Cf. 

THE! 

Cl=< 

C2^R 

C3-( 

F1M 

FLS-- 

FFcr 

F  -•  F  / 
P.EMF 

r  max 

REMF 
RLS= 
RFCC 

WP  n 

ANGt 


F=(#74.g*Z2/(3  .0*PI  )  )*<FEE*FCC-FDD*F3y  )/<FDD*FAf -FCC*FCC  ) 

:off 


E(6,35) 


F<6, '»0)CMfCD,  BE^A,  RFYDU,  CMLS,  CDL  S,CMFF, 
E  =  0.0 
E ( 6, ^5  ) 

00*K=1,MCARD 

*FC*PCE  (K) 

A 1  -  (  (  2  .F*P  I )/  360  )*  AN'  GL  E 

ABS (CCS (THETA1 ) ) ) *C*S ( THETA1 ) 

HC?*(  (UMX**2.0)  /2.0)*D!A*CL 

(P  T*»2  )x:0!  A-:S  lri(  THETA1  )  )  /(  U1X*PEP.) 

CM*C3-( D*CJ ) 

(CMLS*C3-CPl  S*C1  ) 

=(CMFF*C3-CDFF*C1J 

C2 

=ABS (F  )-ABS(Fl ) 

=  Ff,,1M*CN/C2 

=REMI /FM^X 

(AB5 ( F J-ABS (FLS) ) /FMAX  

=  (  ABS(  F)-AL!    (Fl      R))/l  '-'AX 

=TIME+DbLTAl 
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FG 


ST 
EM 
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N'T  IN 
RMAT 
P"AT 
.4,2 
P>  AT 
2X  ,  • 

Cf.«AT 

R.MAT 
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PMAT 

p  '■'  A  T 

LS  «  , 

P'<'AT 

RMAT 

CP 

0 


UE 
UE 
(3 
(' 
Xt 
C 
F  ' 
(  ' 
{■ 

(' 
(• 

12 
(  ' 
(P 


F3.4, 

1 '  t  2  X 
•  UNX= 
0*  ,?X 

,  7X,  • 


0' 
0« 

7X 
0« 

n« 

x, 

0' 
10 


18, F 
»*C! 
•  ,  F  6 
,,T  ! 

FCC  r 
12.4 

,«C<". 
Dl.S  = 
1  2 .  4 

>'    ■    I 

Oh  '  , 


8.4,F 

A=»  ,r 

.4,  2X 

ME/ PC 


,F 
-  i 


9. 


F 
7X 
•C 
6F 
7X 
FF 

6F12.4,4F15 
4) 


10.6,F8.4, T7) 

0.4, 2X,  •  AMP=«  ,  F8.^  ,2X,  »PER  =  «  , 

,  'f;i  =  '  ,Fa.4,2X,  »N."  un=«  ,16) 

R« ,7X, ■ ALPHA1 , 7X, 'CGSA ' , 15X, •SINA 
, • rSINA1 ) 


7X 
•  ,  7  X  , 
,  Zf  12 
ME'  ,9 
12X, 


4  ♦  2  F  1  2  .  4  ) 

, 'C0=« , 7X, • BETA= « ,7X, •  ^EY\0=' , 7X, 

•f MFF=' ,7X, 'COFF^'  ) 

.4  ) 

X, «F« ,9Xf *  Fl*  t9X,J  P.EMF  •  , 9X  ,  ' FL S • , 9X , 

R F  F  R  •  ) 

.4) 
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